Solutions for a massless off-shell two-loop three-point vertex 



On 



Alfredo T. Suzuki and Alexandre G. M. Schmidt 
Instituto de Fisica Teorica - R.Pamplona, 145, Sao Paulo, SP CEP 01405-900, Brazil 

(February 1, 2008) 

Negative dimensional integration method (NDIM) seems to be a very promising technique for 
evaluating massless and/or massive Feynman diagrams. It is unique in the sense that the method 
simultaneously gives solutions in different regions of external momenta. Moreover, it is a technique 
whereby the difficulties associated with performing parametric integrals — the standard approach 
— are transferred to a simpler solving of a system of linear algebraic equations. Employing this 
method, we calculate a massless two-loop three point vertex with all the external legs off-shell. Then 
NDIM approach allows us to obtain twenty-one distinct new power series representations for the 
integral in question. In order to verify the correctness of our results, we consider five particular cases 
where either two of the external legs are put on-shell, or one of them amputaded or one exponent of 
■ the propagators is set to zero, and compare our results thus obtained with the ones calculated with 

' standard methods in positive dimension. 
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ID 1 I. INTRODUCTION. 

^vq I Everyone who has ever studied perturbative quantum field theory (QFT) has come face to face with Feynman 
;> integrals and knows all too well that evaluating these loop integrals is often a hard task. Books on QFT generally 
teach us how to deal with them either in terms of the so-called a-parametrization or by using Feynman parameters for 
the propagators. This approach, contrary to what we may think, is perhaps one of the most laborious techniques to 
solve them because the number of parametric integrals increases with the number of propagators. Moreover, the added 
parametric integrals are often very difficult to solve exactly even if there are no massive particles in the intermediate 
states. There are several techniques that have been developed over the years in order to solve Feynman integrals, 
and we can just mention a few ones: the Mellin-Barnes' representation of massive propagators 0, the Gegenbauer 
polynomial approach in configuration space and integration by parts || and some others [Q. Yet in this arena, the 
simpler the method the better. A novel approach that was suggested some years ago has to do with the use of negative 
JL , dimensions and seems to be a promising technique. 

In mathematical physics, one of the most powerful and useful principles is that of analytic continuation. It is 
the underlying principle that allowed, for instance, 't Hooft et al || to develop the elegant technique of dimensional 
regularization (DREG). Besides a great number of other uses in Physics, there is one that particularly is of interest 
and that most concerns us here: It is that of negative dimensional integration method (NDIM) |6|,[?j] which combines 
the two amazing features in its heart: The principle of analytic continuation and the technique of DREG. What is the 
advantage of such a combination over the plain DREG? The latter one allows us to calculate Feynman diagrams in 
the analytically continued _D-dimensional space, while the former allows us to greatly simplify the technical difficulties 
associated with performing parametric integrals, just by analytically extending D over negative values. This simplifi- 
cation comes about in virtue of the polynomial nature of the integrands. In a few words, this character emerges from 
the equivalence between negative dimensional bosonic integration and positive dimensional fermionic integration |Q. 

What is the price one pays for working them out in negative dimensions? Basically, the technical difficulties that 
arise are solving many systems of linear algebraic equations, performing gaussian/gaussian-type integrals and dealing 
with multi-indexed power series. For the former two, one can just ask whether could it be simpler than these? Yet 
the difficult part is not absent: the laborious piece comes in the form of multi-indexed power series. However, looking 
at it from the encouraging side, we can say that it allows us a standardized representation for Feynman integrals in 
terms of power series. 

In a previous paper || we carried out what we, as theoreticians, call a "lab-test" , that is, with a new approach one 
studies a well-known system. Employing NDIM we calculated a massless two-loop three point vertex, keeping two of 
the external legs on-shell. The NDIM technique led us to discover twelve different ways in which the result could be 
written down, i.e., a twelve- fold degeneracy for that particular integral. We also considered some two-loop self-energy 
diagrams for a massless theory ||. 

Of course NDIM is not the only technique to calculate Feynman integrals. Recently, Fleischer et al studied asymp- 
totic expansions of some two- loop vertex Jlfj] and Frink et al gave results for a general massive two-loop three point 
vertex [||. Ussyukina and Davydychev Jl l| calculated the Feynman diagram we will study in this paper but they 
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did give the result in terms of only two dimensionless variables. On the other hand, we will write down twenty-one 
results in terms of not only two combinations of external momenta, but rather in many different ratio combinations 
for the external momenta . Our aim here is not to make any numerical calculations, of course, but analytical ones. 
The interested reader in the two-loop calculation "technology" in QFT can consult a good review on this subject by 
Davydychev O]. 

The outline for our paper is as follows: in section 2 we calculate the vertex in question in euclidean space, in section 
3 we consider five special cases where only one external leg is off-shell and in the last one, section 4, we conclude the 
work. 

II. OFF-SHELL TWO-LOOP VERTEX. 

This computation is performed following the few simple steps outlined in |^,^). First of all, let us consider the 
integral, 

I = J J d D r d D q exp [-aq 2 - (3(q - p) 2 - 7 r 2 - u(q - r - k) 2 ] , (1) 

which corresponds to the diagram of figure 1. 

The general solution for the integral in negative D, defined by 

Jndim = J J d D q d D r {q 2 f [(q - p) 2 ] 3 (r 2 ) 1 [(r - q + kf] m , (2) 

is given by the multiple series, 

bNDiM = y(i,3,l,m;D) > — j — j — j — j — j — j — j — j — -, (3) 
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where 



g(i, j, i, m; d) = (-7r) D r(i + i)r(i + j)r(i + z)r(i + m)r(i -a- ^d), 

and for convenience we use the definition a = i + j + I + m + D. The system one must solve is, 

ni + n 2 + n 3 + n 5 + n 6 = i 
m + n 2 + 7i4 + nr + n 8 = j 
ni + ria + TI4 + n 5 + n 7 + n g = I (4) 
n 2 + ri3 + ri4 + n 6 + n 8 + n 9 = m 
ni +n 2 + ri3 + n 4 = a. 

It is an easy matter to see that this system is composed of five equations with nine "unknowns" (the sum indices), 
so that it cannot be solved unless it is done in terms of four arbitrary "unknowns" . These, of course will label the 
four remnant summations, which means that the answer will be in terms of a fourfold summation series. From the 
combinatorics, it is a straightforward matter to see that there are many different ways we can choose those four indices; 
indeed, we can choose Cf = 126 different ways. In other words, what we need to do is to solve 126 different systems. 
Of these, 45 are unsolvable systems, i.e., they are systems whose set solution is empty. There remains therefore 81 
which has non-trivial solutions. Of course, the trivial solutions are of no interest at all. However, the non-trivial 
solutions generate a space of functions with different basis, characterized by their functional variable, according to 
the different possibilities allowed for ratios of external momenta. Each basis is a solution for the pertinent Feynman 
integral, which is connected by analytic continuation to all other basis defined by the other sets of solutions. We 
remind ourselves that a basis that generates a given space can be composed of one or several linearly independent 
functions combined in what is called linear combination. 

Each representation of the Feynman integral will be given by a basis of functions generated by the solutions of the 
systems [jl3|| . Of course, only linearly independent and non-degenerate solutions are relevant to define a basis. 

It can be easily seen that the diagram we are dealing with here is symmetric under the exchange of external momenta 
k 2 *-* t 2 . This symmetry is reflected by the systems we have to solve, and the solutions display this fact. Therefore, 
solutions within this category will be given only once, that is, the ones which can be obtained by symmetry will not 
be written down explicitly. 
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A. Two Variables. 



With the solutions in hands — it is an easy matter to write down a computer program to solve all the systems — 
and the general form of the results (^), we can start to build the power series representations of the Feynman graph. 

We begin our analysis of the solutions for the systems by looking at the simpler ones having two variables, defined 
by ratios of external momenta. Of course, there are in fact four sums but two of them have unity argument, making 
it possible for us to actually sum the pertinent series as we shall shortly see. The variables are, 

(x,y), (z^- 1 ), 

5 

where we define the dimcnsionless ratios 

y = ^ (6) 

O 

P 

z = ¥- 

Note that the second pair of variables above is exactly symmetric to the first one by the interchange k 2 «-> t 2 , which 
means that we will not write this second solution explicitly. Also, each of the three pairs above appears twelve times 
among the total of 81 systems with non-trivial solutions. Each of these is therefore twelve-fold degenerate just like in 
the on-shell case calculated in || . A way of expressing the first solution in positive D is given by 

^ n<tln 6 ln 7 ln 9 \ 
(— I — m — ^D\n 4 + n 9 )(^D + i\n 4 + n-? + n 9 ) 

X — 



(7) 



where 



(1 + j + I + \D\ - n 4 + n 6 - m) 

(-1)" 7 ; 

(1 + i — l\m — HQ + n-i + ?^9)(1 — j — I — m — D\nn — no) 

Pi AC = HI - 3 \d){~3\3 + I + \d){-1\1 +m+~D) 

x (-m| - I - \d)(1 + m + D\j)(a + -D\i - a). (8) 

Using the property (a\b+c) — (a + b\c)(a\b) and the well-known summation formula |I^ , [l5|] of Gauss' hypergeometric 
function 2-F1 with unity argument, 

„ , , ,s r(c)T(c — a — b) 

2F1 (0,6; el = )J \ r(, 9 

1 (c — a)L (c — b) 

we can sum the series in uq and in nj above. We then get, 

s^ c = 7r D ( p 2 ) i +3+? D (k 2 y +m+ ? D pA c ( x ) ng (y) n4 



n4,n9— 



(—1 — m — ^D\n4 + ng)(^D + i\n4 + ng) 
(1-j-l-m- D|n 4 )(l + i + j + \D\n 9 ) 

= tt d (p 2 ) l+: > + i D (k 2 ) l+m+ i D Pf c 



X Fa 



-l-m-^D, \D + i 



1 + i + j + \D, 1-j-l-m-D 



x,y\. (10) 
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The remnant double summation in 714 and ng is by definition the Appel's F 4 hypergeometric function jTjJllI of two 
variables. In the particular case where i=j = l = m = — 1 we can simplify even more this result by using a reduction 
formula E3] which relates the F4 function to the gaussian hypergeometric function 2-F1, 



F 4 [a,0;l + a- (3,(3 



—w 



(l-u){l-wY {l-u)(l-w) 
u(l — w) 



1 -u 



For this special case, 



(11) 



D 



2-±D 



rp-3)r 2 (2-|£>)r 3 (iD-i) 



t(d - 2)r(|D - 4) 

2 F, ( 2-i£),i£>-l;4-£> 



/cp p 2 — fcp 



p 2 k 2 — kp 



(12) 



We note that in order to regularize the divergences we can follow the standard procedure of dimensional rcgular- 
ization jl]]. 

The next solution (third one) also gives double series. Following the same steps we can sum two series of unity 
argument and the remaining two are by definition the Appel's hypergeometric function F4, 



s£ c = TT D ( P 2 yp 2 AC 



X F4 ( — (7, — I — m D) 1 + j — <7, 1 + i — a \x 1 ,z 1 



(13) 



where 



P 2 AC = Hk)(-J»(-i| ~m- ^D){-m\l + m+ -D) 



1 



1 



x (l + m + D\-l- ^D)(a + ^D\ -2a - ^D). 



(14) 



Again, in the particular case where all the exponents of propagators are minus one, this F 4 function reduces to a 
gaussian hypergeometric function, too. Using another reduction formula E4], namely, 



F 4 [a,(3; (3,(3 



{i-u){\- w y 

(1 - u) a {\ - w) a 2 F 1 {a, 1 + a - ft (3 \uw) , 



(15) 



we get, 



(p 2 — kp)kp 



-, 2-D/2 



k 2 t 2 



pAC 
r 2 



x 2 Fi \ 2- X -D. X -D-VA-D 



(k 2 -k P y 



kH 2 



(16) 
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B. Three Variables. 



In a manner similar to the previous results, there are solutions which have three remaining variables, meaning that 
one of the series with unity argument is summed out. There are six sets of these, determined by their variables, each 
appearing four times, i.e., a fourfold degeneracy. Just to keep our accounting straight, 4 x 6 = 24 systems yielding 
solutions with three variables. These, added to the 36 of the previous subsection, gives us 60 from the total of 81 
non-trivial systems. 

The solutions within this category have functional dependencies given by: 

(x, x, y), (z, z, y- 1 ) 

(x, y, y), 0, 2T 1 , 2T 1 ) 

(x-\ x~\ z- 1 ), (z-\ z~\ x^) 

Note that the above triplets are conveniently arranged into pairs connected by the symmetry k 2 <-> i 2 , so that in 
the following, only three of them will be dealt with specifically. 

We list below the triple power series representations provided by NDIM, 

3 AC =f ^ (zr^Or 1 )" 3 (m+|gjni) 
3 „ k n n ^ 2 ^ {l-i-3-\D\n,+n 2 ) 



where 



where 



where 



n lt n2,n 3 =0 K ■' 1 

(I + \D\n 2 ){-i\ni + n 2 + n 3 )(-a\ni + n 2 + n 3 ) 
(l + m + D\n\ + n 2 )(l — i — I — m — D\n 3 ) 



h = ir D (t 2 y(-j\<j)(-l\l + m+ \D){-m\l + m+ l -D) 
x (I + m + D\i + j - I -m - ^D)(a+^D\ -2a- ^D), 



,ac f f (lT 1 ) n7+n8 (*) n9 {rn+\D\n 7 ) 

'4 = U 2^ 



n 7 !n 8 !n 9 ! (1 + i + j + \D\n Q ) 

(I + \D\n & ){a + \D\n 7 + n 8 + n 9 )(i + \D\n 7 + n 8 + n 9 ) 
(l + m + D\n 7 + n 8 )(l + i + I + m + D\n 7 + n 8 ) ; 



U = n D (k 2 rxiz i+ i D (-i\ j \D){-j\j + 1+\d) 

x - i-m- D)(-m\i + m + ^D)(l + m + D\-l- ^D), 



^ n4!n 7 !n 8 ! (1 - j + a\n 7 + n 8 ) 



(17) 



(18) 



(m + \D\n 7 ){-j\ni + n 7 + n 8 )(i + \D\n A + n 7 + n 8 ) 
(a — i — j\n 7 + n 8 )(l + i — a\m) 



h = ir D {k 2 y X i(-i\a){-l\l + m+ ^D)(-m\l + m + ±D) 
x (a + ^D\i - a)(l + m + D\i + 2j - 2a). 
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It is not difficult to note that f$ has only one pole in the particular case when i=j = l = m= — 1 and D — 4. 
What is going on here? Looking at the series we observe that there is a factor (1 + i — a\rii) in the denominator. For 
the special case in question it gives, 



T(l + i-a) 



r(4-D) 



(l + i-a\n 4 ) T(l + i-a + n 4 ) r(4 - D + n 4 ) ' 

so that the only term which is not singular is the first one, n 4 — 0, while the others become divergent for D = 4. 
Then, in fact, we have a double pole as it should be. 



C. Four Variables. 



Lastly, we consider in this subsection solutions of systems of linear algebraic equations leading to four variables, or 
better, fourfold summation with four variables. There are 21 of these solutions, which completes the total of 81 non 
trivial solutions for the systems. Again, the functional variables are given paired with their corresponding symmetries, 
as follows: 



(z , z , x , x 

{y, y, z, y- 1 ] 

(x, y, z, y- 1 ] 
(x, y, X' 1 , z~ v 
(z, z, x , z 
(z, y~\ z-\ z~ v 



{x, x, y, y) 
(y- 1 , y-\ x, y) 



(z, y \ z 1 , x *) 

(^X y Z j X ^ 

(x, y, X' 1 , x- 1 ) 



To get the accounting straight, let us again mention that the first three appear just once and the remaining nine 
appear twice, totalling the needed 21 of this category. 

The next two solutions appear just one time; the first one being given by 



where 



and 



{n,}=0 



^-lJns+nB^-l^T+ns (J + ±D\n 6 + Ug) 

n 5 \n 6 ln 7 ln s l (1 - i + a\n 5 + n e ) 

(m + |-D|n 5 + n 7 )(a + \D\n 5 + n 6 + n 7 + n 8 ) 
(1 -j + a\n 7 + n$) 



h = * 



D 



2+2 \ ° 



~p 2 



z l x 3 - l-m-D) 



x - l-m- D)(-l\l + m+ ^D)(-m\l + m + ^D), 



{n,}=0 



(^ni+rw^-ijTiT+ns (-j\ ni +n 2 +n 7 + n 8 ) 



niln 2 ln 7 ln$l (1 - j + a\n 7 + n 8 ) 

{I + \D\n 2 + n 8 )(m + \D\n x + n 7 ) 



{l-i- j -\D\m+n 2 ) 



(20) 



(21) 



where 
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f 7 = n D (k 2 f y 1 -l-m- D){-l\l + m+ ^D) 

x (-m\l + m+^D)(a+-D\i + j-a), 

The next set of five solutions are such that each one of them appears twice but with exponents of propagators 
interchanged, i.e., / <-> m. This means that the space of functions here is generated by two linearly independent basis 
functions, and the series representation for the Feynman integral will be given by 

Sf c = + S^ 2 \ r=8, 9, 10, 11, 12 

where is obtained from S^ Ci ^ by interchanging of the exponents I <-> m. Then, from 



Js ^ n 2 !n 6 !n 7 !n 9 ! (1 + » + j + iD|ng - n 2 ) 1 j 



where 



,AC(1) _ ,(1) g ( y )» 2 +» 6(?/ -l)»v( z )»o (-l)"r+n B(j + + ^ 

(tr - Z|rz. 7 + ng - n 2 ) 
+ l\n 2 + n 6 — n 7 - n 9 )(l + i + m + \D\n 7 - n 6 - n 2 ) 



= (-n) D (t 2 r (iT 1 ) i+m+ * IJ {z) i+ ^ D {-i\ -j l D )(-j\j + l + ±D) 
x (-l\a)(-m\ - i - ±£>)(<7 + -I -a- ^D), 



we have the general solution 



S£ C = S£ C{1) +S* C{2 \ (23) 



where the second term is obtained from the first by interchanging / «-> m. 

In this solution, it is important to note that there is one pole (in the special case when the exponents of the 
propagators are minus one) of the form T(i — I). To regularize it one must introduce a small correction to one of the 
exponents (a suitable one, of course ), that is, to take i = —1 — 6 and then expand all the factors and the power series 
around 5 — 0. Then, this "apparent singularity" cancels out [ pl] , p7| . 

Another result is given by, 

S AC(l) =f (l) y (zT 1 (XT 2 (vT^y- 1 ) 117 (_l)» 1+ » 2 +» 6+ n 7 

9 9 r i n ni!n 2 !n 6 !n 7 ! (1 + i + m + \D\n 7 - n 2 - n 6 ) 



(I + \B\n 2 + n 6 )(m + \D\n x + n 7 ) 



(24) 



where 



so that we get, 



(1 — i — j— \D\nx + n 2 )(l + j + I + \D\n§ — ni — n 7 ) 



jf = ^(fc 2 ) l +"^/ 2 (t 2 )^ i+ ^/ 2 (-z| - m - \d){-j\ -l-\D) 
x + m+ ^D)(-m\l + m+ 7;D)(a + ^D\ -l-m- D), 



S£ c = sf (1) +S* c(2) . (25) 



We note that in this result appears three gamma functions that diverge in four dimensions. What is the nature of 
this extra singularity? It can be a pinch singularity |18| ]. 
Next we have 
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where 



cAC(i) _ ,(i) y- 

{n,}=0 



(i) >^ (z) n >(z- 1 ) n *(x- 1 ) r *(y- 1 ) n * (-!)"'(/ + ^£>|n 2 + n 8 ) 

(1 - j + a\n 7 + n 8 ) 



7i 2 !7i5!n 7 !n 8 ! 



(to + 5-D|n 5 + n 7 )(z + j + to + D|n 5 + n 7 - n 2 ) 
(1 + j + to + |D|n 5 - n 2 - n 8 ) 



/lO } = 7T D (fc 2 )V (P 2 )-™--* (t 2 ) l+1 > D (-i\i + TO + - TO - ID) 



1 



yielding the solution 



Next we have 



x H|2Z + l -D){-m\ -i-l-D){o+ l -D\-l- l -D), 



Q AC _ qAC(1) , Q AC(2) 
D 10 ~ °10 "T" °1 



'10 



S AC W = ,(1) E 



n 2 +n 3 



where 



{£}=a n 2! n 3!«5!«8!(! + * + 1 + ^D\n s - n 3 - n 5 ) 

(Z + jD|n 2 + n 8 )(-j + cr\n$ - n 3 ) 
(1 - j + l\n 2 + n 8 - n 3 - n 5 )(l + j + to + iD|ra 5 - n 2 - n 8 ) ' 



/<}> = ^(p^+'+iW^+^HI - J - ^D)(-i|a) 
x (_/|j)(- m |_ i7 --lD)( <r +IiJ|/- <7 ) ) 



which yields the general solution 



Q AC _ oAC(l) Q AC(2) 
°11 ~ J ll * °1 



'11 



Finally, we have 



jAC(l) _ ,(i) 



(z)™ 2 (y- 1 )" 3 (z- 1 ) n5+n7 
n 2 !n 3 !n 5 !n 7 ! 



7 12 



/12 

{«,}=o 

(-l)" 3+ " 7 (TO+lD|n 5 + n 7 ) 
(1 + z + to + iD|n 7 - n 2 - n 3 ) 

(-/ + ojn 5 + n 7 - n 2 )(-j + ojn 7 - n 3 ) 
(1 + to + cr + \D\n 5 + n 7 - n 2 - rc 3 ) 



where 



A ( 2 } = (-^) D (t 2 Yz- m -^ D (-i\i + m+lD)(-j\a) 



1 



yielding the general solution 



x (-Z|ct)(-to| - i - ^)(ct + - to - 2cr - D), 



oAC _ Q AC(1) q AC(2) 
o 12 — o 12 -|- o-i 



-"12 



D. Region of Convergence. 



The various power series we obtained with negative dimensional integration approach to solve the Feynman integral 
relative to the chosen two-loop vertex diagram are very similar to hypergeometric series. Of course hypergeometric 
functions of more than two variables are known and they are called Lauricella's functions . But in Appel et 
al studied four Lauricella's functions named Fa, Fb, Fc, and Fp — even though they mention that there are other 
fourteen. So, we do not know what are the regions of convergence of them nor even how they are called. However, 
knowing the region of convergence is not an essential thing here, because our external legs arc off-shell anyway. This 
question becomes meaningful in the case where one of the legs are put on-shell, since only in this regime we can 
atribute a value to the external momenta. 



III. ON-SHELL LIMIT. 



Of course, particular cases of on-shell external legs must be contained in the set of off-shell solutions. To check 
on this, let us take two legs on-shell, namely, let k 2 = t 2 = 0. Not all off-shell solutions S AC are suitable for taking 
this particular limit, because some of them either vanish or become divergent. It is easy to see that such a suitable 
solution is given by eq. (|l3|) , because in this limit only the first term in the F4 series is non-zero while all the others 
vanish, leaving us with 

S? c (k 2 , t 2 = 0) = n D (p 2 y(-i\*){-j\a)(-l\ - m - l -D){-m\l + m + l -D) 

x (I + m + D\ - I - ^D){a + -D\ -2a- -D). (32) 

This result is valid for arbitrary D and (negative) exponents of propagators. In order to confront this result with 
known one we still need to go further in specializing to the case where i = j = l = m = — lto get 

52 T{D-2)T{\D-A) ' (33) 

This is the very result obtained by Hathrell using standard procedures for calculating Feynman diagrams in 
positive D. Of course, a more straightforward way of getting this result using NDIM is to put the corresponding legs 
on-shell from the very beginning, and what we get then is twelve systems to solve with non-trivial solutions, exactly 
the number we have for the solution in question: a twelvefold degeneracy giving the same correct result B . 

We can consider also other special case, namely, the one where p 2 = k 2 = 0. This one is interesting because it 
contributes to other two-loop three-point diagram |2p[| if we apply the integration by parts technique |j| . The general 
result, for arbitrary (negative) exponents of propagators and (positive) dimension can be read from the solution S3, 
eq.©, 

S£ c {p 2 , k 2 = 0) = ir D (t 2 ) a (-j\a){-l\l + m+ ^D)(-m\l + m + ^D) (34) 

x (l + m + D\i+j -l-m - ~D)(<r + -D\ -2a- ^D), 
taking the same particular case of Kramer [ pp[ et al, i.e., I = m = —2 and i = j = —1, we obtain, 

sic = ^ ^-^-5)r(6- jD)> 



which is the well-known result in euclidean space. 

Two other simpler special cases can be read from this graph, namely, when p — 0, k — t (see fig. 2) and k = 0, p = t 
(see fig. 3). The solution S3 gives us the first one, 



S^(j) = 0,k = t) = f 3 (-i-j\i)(i + l + m + D\-i), (36) 

note that in the n\ and ni series only the first term contributes and the n.3 one reduces to a gaussian hypergeometric 
function with unit argument. 







The second case (see fig. 3), k = 0, p = t, can be read from 52, 

S£ c (k =0,p = -t) = K D (p 2 yP 2 AC (-i +a\- *){\D+j - a^), (37) 

this result can be used to study two self-energy two-loop graphs and agrees with our previous results || . 

Finally, let us check up on a solution that has four series. Let j = 0, we get the graph of figure 4. The solution 
that allows us to consider this limit is SV, 

Sf c {j = 0) = n D (k 2 rf r (j = 0), (38) 

observe that there is no sum in the result. This is due to the factor (—j\ni + n<i + n-j + n&) which leads to only one 
non- vanishing term, i.e., when n\ — = n-i = n§ = in the series. 



IV. CONCLUSION. 



We have shown in this paper how we can work out a two-loop vertex diagram with all external legs off-shell using 
the NDIM technique to solve the pertinent Feynman integral. Altogether, twenty-one distinct results are obtained 
via NDIM technique for the considered two-loop three-point vertex diagram. These are expressed in terms of power 
series which can be identified as hypergeometric functions. The simpler ones are Appel's F4 hypergeometric function 
with two variables, which for the particular cases where all the exponents of the propagators are set to minus one, 
can be transformed into even simpler ones of the gaussian hypergeometric function type. The technique provides 
simultaneously with several analytic continuation formulas between different results, because they arise from the 
same Feynman integral (Q). 
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Figure 1: Off-shell two-loop three point vertex calculated with NDIM. 



1 




Figure 2: A special case ( p = 0, k = t ) of the diagram of figure 1. 



Figure 3: Taking k = and p = t the graph of figure 1 reduces to this one[9]. 
The former we get setting i=j = l = m = — lin eq.(37) and the latter 
making i = — 2 and the others equal to minus one in the same equation. 



1 




2 




Figure 4: A special case ( j 



= ) of the diagram of figure 1. 



